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A comprehensive analytical model of the interaction of the cantilever tip of the 
atomic force microscope (AFM) with the sample surface is developed that accounts 
for the nonlinearity of the tip-surface interaction force. The interaction is modeled as 
a nonlinear spring coupled at opposite ends to linear springs representing cantilever 
and sample surface oscillators. The model leads to a pair of coupled nonlinear 
differential equations that are solved analytically using a standard iteration procedure. 
Solutions are obtained for the phase and amplitude signals generated by various 
acoustic-atomic force microscope (A-AFM) techniques including force modulation 
microscopy, atomic force acoustic microscopy, ultrasonic force microscopy, 
heterodyne force microscopy, resonant difference-frequency atomic force ultrasonic 
microscopy (RDF-AFUM), and the commonly used intermittent contact mode 
(TappingMode) generally available on AFMs. The solutions are used to obtain a 
quantitative measure of image contrast resulting from variations in the Young 
modulus of the sample for the amplitude and phase images generated by the A-AFM 
techniques. Application of the model to RDF-AFUM and intermittent soft contact 
phase images of LaRC-cp2 polyimide polymer is discussed. The model predicts 
variations in the Young modulus of the material of 24 percent from the RDF-AFUM 
image and 1 8 percent from the intermittent soft contact image. Both predictions are in 
good agreement with the literature value of 21 percent obtained from independent, 
macroscopic measurements of sheet polymer material. 


PACS numbers: 68.37.Tj, 81.07.-b, 82.35.Np, 68.37.-d 
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I. INTRODUCTION 

The atomic force microscope 1 (AFM) has become an important nanoscale 
characterization tool for the development of novel materials and devices. Dynamic 
implementations of the AFM (we shall call acoustic-atomic force microscopies or A- 
AFM)) such as intermittent contact mode (TappingMode), force modulation microscopy 2 
(FMM), atomic force acoustic microscopy 3 ' 4 (AFAM), ultrasonic force microscopy 5 ' 6 
(UFM), heterodyne force microscopy ' (HFM), resonant difference-frequency atomic 
force ultrasonic microscopy 9 (RDF-AFUM) and variations of these techniques 10 ' 14 utilize 
the interaction force between the cantilever tip and the sample surface to extract 
information about sample material properties. Such properties include sample elastic 
moduli, adhesion, surface viscoelasticity, embedded particle distributions, and 
topography. The cantilever tip-sample surface interaction force is generally nonlinear 5 , 
although in some operational modes the interaction force can be taken to a good 
approximation to be linear. A comprehensive treatment of the interaction force, however, 
is lacking because of the difficulty in accounting for the nonlinear terms. We consider 
here a detailed analytical treatment of the cantilever tip-sample surface interaction force 
that includes the lowest-order terms in the nonlinearity. Such terms are sufficient to 
account for all operational characteristics and material properties obtained from the 
various acoustic-atomic force microscopies cited above. 

We begin in Section II by developing a realistic mathematical model of the 
interaction between the cantilever tip and the sample surface that involves a coupling via 
the nonlinear interaction force of separate dynamical equations for the cantilever and the 
sample surface. A general solution is found that contains static terms (including static 
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terms generated by the nonlinearity), linear oscillatory terms, and nonlinear oscillatory 
terms. Individual or various combinations of these terms are found to apply as 
appropriate to a description of a particular acoustic-atomic force microscopy. The 
complexity of the general solution begs consideration of ways to simplify the 
mathematical expression by estimating the relative magnitudes of the material and 
dynamical parameters embedded in the equation. Such estimates are considered in 
Section III. Application of the solution to each of the above-cited acoustic-atomic force 
microscopies is given in Section IV using simplifications resulting from the estimate of 
parameter values obtained in Section III. Section V provides an analytical analysis of 
image contrast for each of the A-AFM techniques addressed in the model. Application of 
the model to RDF-AFUM and intermittent soft contact phase images of LaRC-cp2 
polyimide polymer is discussed in Section VI. 

II. ANALYTICAL MODEL OF NONLINEAR CANTILEVER DYNAMICS 
A. General dynamical equations 

The cantilever of the AFM is able to vibrate in a number of different modes in 
free space corresponding to various displacement types (flexural, longitudinal, shear, 
etc.), resonant frequencies, and effective stiffness constants. Each cantilever mode n of a 
given displacement type may be represented by an effective mass nic attached to a spring 
of stiffness constant k cn and resonant frequency a> cn = ^k cn / m c . Since each mode is 

subject to the same driving forces, we may express the cantilever time-dependent 
displacement tj c in its most general form as the sum over all modes and write 
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7c(0 = E7c«(0 (1) 

n 

where r] cn is the cantilever displacement corresponding to mode n and t is time. During 

contact of the cantilever with the sample surface, vibrations of the cantilever give rise to 
oscillations of the sample surface via the cantilever tip-sample surface interaction forces. 
Likewise, an oscillating sample surface resulting, for example, from an incident 
ultrasonic wave generated at the opposite surface of the sample will also give rise to 
oscillations of the cantilever via the same interaction forces. The vibrating sample 
surface may be represented by an effective mass m s attached to a spring of stiffness 

constant k s . 

For definiteness we consider only flexural modes of the cantilever and normal 
out-of-plane oscillations for the sample surface as indicated in Fig.l. The displacement 
of the cantilever tip is tj c and the normal displacement of the sample surface is77 s . The 

quiescent, equilibrium distance between the cantilever tip and the sample surface in the 
presence of an interaction force is zo. The tip-surface separation distance at an arbitrary 
time, when the cantilever or sample surface or both are in oscillation in the presence of 
the interaction force, is z. The spring model representing the dynamics of the tip-surface 
interaction is also shown in Fig.l. Both the cantilever and sample surface springs are 
assumed to be linear with stiffness constants k cn and k s , respectively, while the 

nonlinear interaction forces are represented by a nonlinear spring having a linear stiffness 
constant F' and a nonlinear stiffness constant F" as indicated in the figure. 
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When the cantilever tip is in contact with the sample surface, the nonlinear tip- 
surface interaction force F(z) provides a coupling of the cantilever and surface 
oscillations. The dynamics may be represented by a set of coupled differential equations 
for each mode n as 


m cVcn + Ycicn + k cnVcn = F ( z ) + F c COS (O c t (2) 

™ s V sn + YsVsn + KVsn = F ( z ) + F s cos (p) s t + 0) (3) 

where rj cn (positive down) is the cantilever tip displacement for mode n, rj sn (positive up) 
is the sample surface displacement for mode n, co c is the angular frequency of the 
cantilever oscillations, cu s is the angular frequency of the sample surface vibrations, y c is 
the damping coefficient for the cantilever, y s is the damping coefficient for the sample 
surface, F c is the magnitude of the cantilever driving force, F s is the magnitude of the 

sample driving force that we assume here to result from an incident ultrasonic wave 
generated at the opposite surface of the sample, and 6 is a phase contribution resulting 
from the propagation of the ultrasonic wave through the sample material. 

We have shown previously 9 that for an acoustic wave propagating through a 
sample of thickness a/2 with phase velocity c and wave number k, containing an 
embedded feature of thickness d/2 for which the phase velocity is ca, the total phase 
contribution 0 is given by 
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where 


0 = -(% + Az) 


ka _i sin ka 

Z = — + tan 

2 e aa - cos ka 


( 4 ) 

( 5 ) 


A Z = -¥ 


1 

— + ■ 


e aa cos ka - 1 


] 


(e aa - cos kay + sin z ka\j 


( 6 ) 


and 


¥ = m {C(i c) . (7) 

c d 

The factor - j is the contribution to the phase from the featureless bulk material and 
- A z is the contribution from a phase variation due to the embedded feature. 

A typical nonlinear interaction force F(z) is shown schematically in Fig.2. This 
force results from a number of possible fundamental mechanisms including electrostatic 
forces, van der Waals forces, interatomic repulsive (Bom-Mayer) potentials, and Casimir 
forces 15 . It is also influenced by chemical potentials as well as hydroxyl bonds resulting 
from atmospheric moisture accumulation on the cantilever tip and sample surface 16 . 

We note from Fig.l that for a given mode n, z = z a - ( rj cn + r/ sn ) . We use this 

relationship in a power series expansion of F(z ) about z Q to obtain 


F(z) = F(z 0 ) + F' (z 0 )(z - z 0 ) + ^ F” (z 0 )(z - z 0 ) 2 + • 


( 8 ) 
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F(. z q) F (z Q )(?}cn 2 ^ ( z 0 X^cn ) "*" ' 


where the superscripted prime denotes derivative with respect to z. Substitution of 
Eq.(8) into Eqs.(2) and (3) gives 


m cVcn YcVcn [^c/i ^ (• z o)]^7cn (.^o'PJsn -^'(• z o) COS<W c i 


( 9 ) 


■*" 2 ^ (^0 )(n cn ^1 sn) ' 


™sVsn Y sV sn ^ ( z o)^7sh ^ ( z 0 )^7cn ^X^o) cos(<w s t + 0) 


+ 2 ^ ( z 0 X*7 cw + ) + ' 


( 10 ) 


Equations (9) and (10) are the coupled equations representing the cantilever tip-sample 
surface dynamics resulting from the nonlinear interaction force at the sample surface. 

B. Solution to general dynamical equations 

We solve the coupled nonlinear Eqs.(9) and (10) for the steady-state solution by 
writing the coupled equations in matrix form and using a common iteration procedure to 
solve the matrix expression. The first iteration involves solving the equations for which 
the nonlinear terms are neglected. The second iteration is obtained by substituting the 
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first iterative solution into the nonlinear terms of Eqs.(9) and (10) and solving the 
resulting equations. The procedure provides solutions both for the cantilever tip and the 
sample surface displacements. Since the procedure is much too lengthy to reproduce here 
in full detail, only the salient features of the procedure leading to the steady state solution 
for the cantilever displacement rj c = £ r/ cn are given. We begin by writing 


Vcn 


cn 


and 


Vsn £ sn 


sn 


( 11 ) 

( 12 ) 


where e cn and % cn represent the first iteration (i.e. linear) static and oscillatory solutions, 
respectively, for the nth mode cantilever displacement, £ cn represents the second iteration 
(i.e., nonlinear) solution for the nth mode cantilever displacement, and e sn ,^ sn , and £ sn 

are the corresponding first and second iteration nth mode displacements for the sample 
surface. 

i. First iterative solution 

The first iterative solution is obtained by linearizing Eqs.(9) and (10), writing the 
resulting expression in matrix form, and solving the matrix expression assuming 

sinusoidal driving terms F c e ia>ct and F s e ia>st for the cantilever and sample surface, 

respectively. The first iteration yields a static solution e cn and an oscillatory solution 


% cn for the cantilever. The static solution is given by 
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k s F ( z o ) 

k cn k s + F '( z o)( k cn +k s ) 


( 13 ) 


The first iterative oscillatory solution is given by 


*3 c« Qcc cos(<y c t + & cc $cc ) Qcs cos (o) s t (f> ss + 0) (1 4) 

where 

<!>cc = ( 15 ) 

tan -l G>c(.Ys k cn + Yc k s)-<°l(rs m c + Yc m s) +F \ z o) a >A.Ys + Y C ) 

k cn k s + m s m c G) c ~ Q^s k cn ™c k s Y cY s) F ( z o)( k cn k s ~ ~ ^c®c ) 


</>ss= ( 16 ) 

tan -l ( Ys k cn + Yc k s ) - a>l ( 7s m c + Yc m s ) + F '( z o ) a s (Ys + fc ) 

k cn k s — (m s k cn + m c k s + y c y s ) + ^ ( z o)(^cm + k s ~ ~ ) 


Qcc = F ci\- k s ■*" F ( z o) ~ m s®c ] Ys a c } {[ k cn k s m s m c®c 

— co c (m s k cn + m c k s + y c y s ) + F (z 0 ){k cn + k s — m s co c — m c co c )] ( 1 7 ) 

[® c (rAcn Yc k s ) — ®c (Ys m c Yc^s ) F ( z o^®c(Ys 7c)] } 
and 

fics ~ F s F (. z o){\- k cn k s + Ms^c^s ~ ( m s k cn +m c k s + YcYs^ 
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+ F ( z 0 )(k cn + k s Ms^s m c^s )] ( 18 ) 

\®s (y s^cn Yc^s ) — (.Y s™c Yc^s^ "*" F ( z o (.Ys Yc )] } 

ii. Second, iterative solution 

The second iterative solution £ cn for each mode n of the cantilever is 

considerably more complicated, since it contains not only sum-frequency, difference- 
frequency, and harmonic-frequency components, but linear and static components as 
well. The second iterative solution £ cn is thus written as 


Ccn G cn,stat + G cn,lin + C cn,diff + C cn,sum + C cn,harm 0 9) 

where Ccn,stat I s a static or “dc” contribution generated by the nonlinear tip-surface 
interaction, C cn ,iin I s a generated linear oscillatory contribution, C C n,diff I s a generated 

difference-frequency contribution resulting from the nonlinear mixing of the cantilever 
and sample oscillations, C cn ,sum is a generated sum-frequency contribution resulting from 

the nonlinear mixing of the cantilever and sample oscillations, and C cn ,harm are 
generated harmonic contributions. 

Generally, the cantilever responds with decreasing displacement amplitudes as the 
drive frequency is increased above the fundamental resonance, even when driven at 
higher modal frequencies. Thus, acoustic-atomic force microscopy methods do not 
generally utilize harmonic or sum-frequency signals. For expediency, such signals from 
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the second iteration will not be considered here. Only the static, linear, and difference- 
frequency terms from the second iteration solution are relevant to currently-used A-AFM 
modalities. 

The static contribution generated by the nonlinear interaction force is obtained to 
be 




cn,stat 


1 k s F"(z 0 ) 

4k C nk s +F'(z 0 )(kcn+k s ) 


[2 £ o + Qcc + Qcs + Qsc + Qss 


+ 2j QccQsc C0S ( a cc ^4*cc ) 2 QcsQss COS GC ss ] 


( 20 ) 


where 


i^cn k s ) F(z 0 ) 

^cri ■*" F ( z o )(k C n + k s ) 


( 21 ) 


Qsc F C F ( z o){\ k cn k s + m s m c co c co c (tn s k cn + tn c k s + y c y s ) 

+ F'{z 0 )(k cn +k s - m s (Qc ~ m c C0 l )] 2 (22) 

[®c (y s^cn y c^s ) — {.Ys^c Yc^s ) F ( z o )®c (Ys Yc )] } > 

Qss = Fs {[kcn F (z 0 ) — m c co s ] + y c co s } {[k cn k s + m s m c co s 

— G) s ( m s k cn + wicks Y cY s') F {z 0 )(k cn + k s — m s 6) s —m c (O s )] (23) 

[®s (.Yskcn Ycks) ~ ®s (.Ys^c Yc^s ) F (z Q )(Q S (y s + y c )] } , 
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a cc = tan 


-1 Ys°>c 

k s +F'(z 0 )-m s G) t 


2 


(24) 


a ss = tan 


-1 Yc m s 


k C n +F \z 0 )-m c a) t 


2 ’ 


(25) 


and <f) cc is given by Eq.(15), Q cc by Eq.(17) and Q cs by Eq.(l 8). 

The linear oscillatory contribution Ccn,lin generated by the nonlinear interaction 
force in the second iteration is obtained to be 


Ccn,lin D ^cjF' ( z o)[ Qcc Qsc ^QccQsc COSflT cc I COS {(O c t '2</> cc + P c + fi cc ) 

'Vc 

(26) 

+ — EqF (z a )\ Q ss + Q cs + 2Q SS Q CS coscc ss ] cos (co s t — 2</> ss + P s + fi ss + 0) 

^ss 

where 


Mcc = tan 


Qcc sma cc 
Qcc cosa CC + Qsc 


(27) 


Mss = ten 


Q ss sin a 55 
Qss cosa SS + Qcs 


(28) 


P c = tan 


Y s <°c 


k s rn s a> c 


(29) 
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A = tan" 1 W , (30) 

^5 — m S^S 

D c =[(k s -m s <DZ) 2 +rW?' 2 . (31) 

A=[(* s -'»A 2 ) 2 +r 2 ® 2 ] 1 ' 2 . (32) 

R ss — ( 1 + ytigtnpOg 0) s ^fn s k crl + tn c k s ^ YcY s) ^ ^ (^o)(^cn ^ k s ~ m s 6) s tn ( jj) s )] 

\Y®s)Ys^cn "*" Yc^s) ~ ®s (Y s***c "*" Yc^s) F ( z o)®s(Ys "*" Yc)) i (33) 

and 

Rcc ~~ { \- k cn k s m s m c<°c — °>c ( m s k cn ™c k s YcYs ') Jr ^ C^o )( k cn k s ~~ ™s®c ~ ^c®c )] 

+ [<»c(Ys k cn + Y c k s ) - <°l ( r s m c + Y c m s) + ^^o^ciYs + Yc)f} 1 ' 2 ■ ( 34 ) 

The difference-frequency contribution Ccn,diff generated by the nonlinear 
interaction force in the second iteration is obtained to be 

C cn,diff — eos[(<y c — co s )t — <j> cc + tf) ss + (3 CS — </> cs + r — $] (35) 


where 



15 


G n = --f- F"(z 0 ){QlQl + QlQl + Qc C Qss + QcsQsc (36) 

^ ^CS 

2Q C c Qcs Qsc Qss COS (pc cc + &ss ) ^QccQcsQss COS OT ^ 

IQccQcsQsc cosa cc ^QscQssQcs cosa ss 
20 cc QssQcsQsC C0s ( a cc — & SS )) ’ 



^csl k cn k s m s k cn(®c ) ^c^s(®c (39) 

— YcYs^c ~ ®j) ^ ( z o)[^cn — *^s(®c — ®s) ~ m c(.®c~®s) ] > 


^C,S'2 = (®c - X/A + Yc k s)- (°>c - <°s ) 3 ) (40) 

+ F'(z 0 )(<y c -® s )(/ s + / c ) , 



( 42 ) 
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and 


F — tan -1 QccQcs sinof cc QscQss^^^ss QccQss^^P^cc ^ ss ) 

QccQcS COS & cc QscQsS COSC?^ + QccQsS C °s(^cc — ^ ss ) QcsQsc 


. ( 43 ) 


Hi. Salient features of the solution set 

The total static solution to the coupled nonlinear equations (9) and (10) for the 
cantilever i) cn stat is the sum of the contribution s cn , given by Eq.(13), from the first 

iterative solution and the contribution f cn ^ tat , given by Eq.(20), from the second 
iteration as 


Vcn,stat £ cn + Ccn, 


stat 


(44) 


The total linear solution T] cn Un to Eqs.(9) and (10) is the sum of the contribution <^ cn 
given by Eq.(14) and the contribution C cn ,li n given by Eq.(26) as 


Vcn.lin £cn + Ccn,lin 


(45) 


The total difference-frequency solution ricn,diff to Eqs.(9) and (10) is simply the 
contribution Ccn,diff S iven b Y Eq,(35). 

It is interesting to note that e cn and the e Q component in ij cn stat do not 
explicitly involve the cantilever drive amplitude F c and the sample surface drive 


amplitude F s , although other terms involving the Q factors, given by Eqs.(17), (18), (22), 
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and (23), in C cn ,stat do involve these drive amplitudes. This means that only the 

contributions stemming from the nonlinearity in the cantilever tip-sample surface 
interaction force respond directly to variations in the drive amplitudes and in particular to 
the physical features of the material giving rise to variations in F s . Further, the 

magnitude of all second iteration (i.e. nonlinear) contributions, Ccn,stat , Ccn,lin’ anc * 
Ccn,diff are strongly dependent on the cantilever tip-sample surface separation z 0 , since 
the value of the nonlinear stiffness constant F"(z 0 ) that dominates these contributions is 
highly sensitive to z 0 . Indeed, F"(z 0 ) attains a maximum value near the bottom of the 
force-separation curve of Fig.2. 

It is interesting to note that the condition k cn k s + F'(z 0 )(k cn + k s ) = 0 gives a 
singularity in Eq.(13) for e cn , in Eq.(20) for Ccn,stat > an d in Eq.(26) [via e 0 of Eq.(21)] 
for C cn ji n . This condition on F'{z a ) does not admit a solution to the original coupled 

Eqs.(9) and (10), since such a condition leads to a value of zero for the secular 
determinant of the coupled equations. On the other hand, the conditions F'(z 0 ) = -k cn 

or F'(z 0 ) = -k s do give rise to valid solutions to Eqs.(9) and (10) for all driving forces, 
static and non-static. The negative value of F'(z 0 ) in this case means that the tip-surface 
separation distance z 0 is smaller than the separation corresponding to that at the bottom 
(absolute minimum) of the force-separation curve. Eqs.(9) and (10) show that a negative 
value of F'(z 0 ) reduces the effective magnitude of the stiffness constants both for the 

cantilever and the surface displacements. Indeed, the values F'(z 0 ) = -k cn and 
F'(z 0 ) = -k s give rise to effective null stiffness constants for the cantilever and surface 
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displacements, respectively, as indicated by the coefficients of i) cn and rj sn in Eqs.(9) 
and (10). 

A further decrease in tip-surface separation z a produces even larger negative 
values ofF'(z 0 ). From Eq.(9) we may define the effective cantilever resonance 
frequency for mode n as (D cn , e ff = silken + FX z o)V m c . Thus, when F'(z 0 )<-k cn , 
o) cn ,eff becomes purely imaginary and produces an additional damping factor in the 

solutions to Eqs.(9) and (10). Such damping is reported by Rabe, Janser, and Arnold 17 . 
It is straightforward to obtain mathematically the effects of such an additional damping 
factor for the cantilever displacement amplitudes Tj cn , if one considers only Eq.(9) and 

ignores the coupling term rj sn . For a driving force given by F c e ic ° ct the magnitude of the 

steady state solution to the equation is | r] cn |= F c {[k cn + F'(z {) ) - m c a)^] +y c co^} . 
For a fixed drive frequency co c , effective cantilever mass m c , and cantilever spring 
constant k cn , the factor F'(zq) produces a decrease in the cantilever displacement 

amplitude | rj cn | when | F'(zq) |>| k cn -m c co c |. Note that it matters little whether 
F'(zq ) is positive or negative (i.e., whether z 0 is greater than or less than the value of z 0 at 
the bottom of the force-separation curve). As the cantilever is brought ever closer to the 
sample surface F'(zq) becomes ever larger in magnitude and \i) cn \ for each mode n 

continues to decrease until the repulsive force in the F(z a ) curve of Fig.2 exceeds the 
fracture strength of the cantilever. 

Finally, it is important to note that for large deflections of the cantilever that 
generally occur for hard contact (linear regime), large bending moments are introduced 
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that produce significant frequency shifts in the cantilever resonance frequencies quite 
apart from those introduced by the interaction force stiffness constant F'(z 0 ) ■ For the 

assessment of F'(z 0 ) near the bottom of the force-separation curve where the 
nonlinearity F"(z 0 ) is maximum (nonlinear regime) and F'(z 0 ) is relatively small, a 
reasonable estimate of F'(z 0 ) can be obtained directly from differences in the contact 
and non-contact resonance frequencies of the cantilever. For the assessment of F'(zq) in 

hard contact (linear regime) it is necessary to account for the bending moments of the 
cantilever in large deflection. 

III. ESTIMATES OF PARAMETER MAGNITUDES 

The results obtained in Section II are general equations of the cantilever response 
resulting from the cantilever tip-sample surface nonlinear force interaction. The 
equations are valid for all materials and material systems, including bio-materials. For 
many materials and microscope operating conditions, the general equations can be 
simplified by using estimates of the relative sizes of the parameters appearing in the 
equations. It is especially useful to obtain estimates of the order of magnitude of the 
parameters m c , m s , y c , y s , k cn , and k s that appear in the model. The estimates and 

resulting equation simplification are, of course, highly dependent on the particular 
material under investigation and on the specifications of the cantilever itself. 

The cantilever stiffness constant k c \ corresponding to the fundamental mode (n = 
1) oscillation is generally measured to lie in the range 0.02 - 80 N m' 1 . A calculation of 
m c can be obtained from measurements of the fundamental resonance angular frequency 
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2 

of the cantilever in free space co c \ and k c \ from the expression m c = k c \ / a> c \ . 
Measurements 17 of k c \ and co c \ for typical commercial cantilevers indicate that m c is of 
the order 10' 11 - 10' 12 kg. 

We may obtain a crude estimate of m s by noting that an oscillating cantilever tip 

in “point contact” with the sample surface generates spherical waves of frequency co 
with pressure amplitude 


p _ P o c i(a>t-k r r) 
r 


(46) 


where P 0 is the pressure at unit radius, r is the radial distance from the contact point, and 
k r is the radial wave number. We assume for definiteness that the “unit” radius is equal 
to the wave oscillation amplitude at the point of contact. Since “point contact” at the 
atomic level ideally involves an interaction between a single atom at the cantilever tip 
and a single atom of the sample surface, we assume that the wave displacement at the 
contact point corresponds to the relative motion between the two atoms. Assuming that 
the relative motion is of the order of interatomic distances, we estimate that the wave 
displacement amplitude is roughly 0.1 nm. From Eq.(46) we assume that when r = 100 
“unit” radii, the magnitude of p is sufficiently small that the oscillating mass associated 
with such small pressures can be neglected. For a unit radius of 0.1 nm we obtain r = 
10' 8 m. It is of interest to note that this value of r is in agreement with estimates 14 of the 
cantilever-sample contact radius r c obtained from Hertzian contact theory using a contact 
force of 200 nN and a cantilever tip radius of 100 nm for a typical polymer material. The 



21 


hemispherical volume V of material corresponding to this radius is of the order 10' 24 m 3 . 
Assuming a mass density p of order 10 3 kg m' 3 for the sample material (roughly a typical 
polymer density), we obtain m s = pV ~ 10' 21 kg. 

The magnitude of the sample stiffness constant k s and the sample absorption 
coefficient y s can be estimated from the relationships of these parameters to the 
longitudinal elastic modulus and acoustic attenuation coefficient of the material, 
respectively. Consider a periodic solid lattice with lattice spacing L between lattice 
points of mass m. Let u n represent the displacement of the nth lattice point from its 
equilibrium position resulting from a force F n acting at that point. We assume that F n 
is a function both of u i and the particle velocity u t where i takes the values n, n-1, and 

n+1. Hence, the force depends not only on the displacement and particle velocity of the 
nth lattice point but also on the displacements and particle velocities of points to either 
side of the nth point. Following the approach of reference 18 we may write lfom 
Newton’s Law 


m 


d 2 u 
dt 7 


n _ 


F n = k s K«n+ 1 “ “n ) “ (“n " “n-1 )] + 7s [(“n+1 “ “n ) “ (“n “ “n-1 )] ( 47 ) 


where k s is the sample stiffness constant and y s is the sample or lattice absorption 
coefficient. 

We divide Eq.(47) by AL = V, where A is some cross sectional area of the lattice 
and V is the volume enclosed by A and L. Taking the continuum limit as L approaches 
zero, we obtain the wave equation 
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(P u 2tPu 

_ 2 ~ C ~ 2 + ^dis 
dt L dx, L 


ePu 

dtdxP 


(48) 


1 12 

where c = (C\\ Ip) is the sound velocity, p = m/AL is the mass density of the material, 
Q i=k s LI A is the longitudinal elastic modulus, an &h dis = y s LI pA is the damping or 
dissipation coefficient for the continuous wave. The magnitude of k s may be estimated 
from the relation k s = C\\A/L where A is the surface area subtended by the above- 
considered hemisphere of radius r. Assuming C\ j * 10 9 -10 n Pa, r « 10' 8 m, and L « 0.4 
nm, we obtain k s » (10 2 - 10 4 ) N m' 1 . 

We may estimate the magnitude of y s by assuming a solution to Eq.(48) of the 

form u = UQe 1 ^-®^ where k is a complex wave number. Substituting this form into 

2 2 2 —1 

Eq.(48), we obtain k =(co/c) [1 - i(o)h dis !c ] . The imaginary part of k is the 
acoustic attenuation coefficient a from which we obtain 


Ys = 


ISpPa 
La 2 


(49) 


An examination of the attenuation coefficients of a variety of materials in both the liquid 
and solid states together with other parameters occurring in Eq.(49) reveals that the 
magnitude of y s for most materials falls in the range 10' 4 - 10' 7 kg s' 1 . 
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IV. APPLICATION OF MODEL TO VARIOUS ACOUSTIC-ATOMIC FORCE 
MICROSCOPE TECHNIQUES 

The estimate of parameters obtained in Section III may be used to simplify the 
equations derived in Section II describing the cantilever response resulting from the 
interaction with the nonlinear cantilever tip-sample surface forces. The specific 
simplification depends on the particular acoustic-atomic force microscopy under 
consideration, the frequencies employed, and the material under investigation. We shall 
consider the most frequently used A-AFM modalities including resonant difference- 
frequency atomic force ultrasonic microscopy (RDF-AFUM), heterodyne force 
microscopy (HFM), ultrasonic force microscopy (UFM), atomic force acoustic 
microscopy (AFAM), force modulation microscopy (FMM), and the most commonly 
used intermittent contact mode. 

A. Resonant difference-frequency atomic force ultrasonic microscopy and heterodyne 
force microscopy 

Resonant difference-frequency atomic force ultrasonic microscopy (RDF-AFUM) 
employs an ultrasonic wave launched from the bottom of a sample, while the cantilever 
of an atomic force microscope, driven at a frequency differing from the ultrasonic 
frequency by one of the contact resonance frequencies of the cantilever, engages the 
sample top surface. It is important to note that at high drive amplitudes of the ultrasonic 
wave or cantilever (or both) the contact resonance frequency generating the difference- 
frequency signal may correspond to one of the nonlinear oscillation modes of the 
cantilever. As pointed out in Section II.B.hi, the effective contact cantilever resonance 
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frequency for the (linear or nonlinear) mode n, ®cn,eff > ma Y be defined as 
oi>cn,eff = s]\ k cn + F*(z 0 )]/m c > where k cn is the cantilever stiffness constant 
corresponding to the nth (linear or nonlinear) non-contact resonance mode. Since F\zq) 
is negative at the separation distance zo corresponding to the maximum nonlinearity, i.e. 
maximum F"(z 0 ) , the resonance frequency of the cantilever in contact is smaller than the 

value when not in contact. The nonlinear mixing of the oscillating cantilever and the 
ultrasonic wave in the region defined by the cantilever tip-sample surface interaction 
force generates difference-frequency oscillations at the cantilever contact (linear or 
nonlinear) resonance. 

Variations in the amplitude and phase of the bulk wave due to the presence of 
subsurface nano/microstructures as well as variations in near-surface material parameters 
affect the amplitude and phase of the difference-frequency signal. These variations are 
used to create spatial mappings generated by subsurface and near-surface structures. 
Heterodyne force microscopy (HFM) also utilizes difference-frequency signals generated 
by the nonlinear mixing in the cantilever tip-sample surface interaction region. In this 
technique no special advantage is taken of cantilever resonances and the difference- 
frequency utilized is generally well below that of the cantilever resonance. 

In both RDF-AFUM and HFM the cantilever difference-frequency response is 
obtained from the nonlinear mixing in the region defined by the tip- surface interaction 
force. The interaction force varies nonlinearly with the tip-surface separation distance. 
The deflection of the cantilever obtained in calibration plots is related to this force; for 
small slopes of the deflection versus separation distance, the interaction force and 
cantilever deflection curves are approximately related via a constant of proportionality. 
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The maximum difference-frequency signal amplitude occurs when the quiescent 
deflection of the cantilever approaches the bottom of the force well, where the maximum 
change in the slope of the force versus separation curve (hence maximum interaction 
force nonlinearity) occurs. 

The dominant term or terms for the cantilever difference-frequency displacement 
in Eqs.(l) and (1 1) depend on the values of k cn both for the linear and the nonlinear non- 

contact modes of cantilever oscillation, A a> , and the value of F'(z 0 ) obtained at the 

separation distance zo at which the maximum difference-frequency signal occurs. We 
designate the non-contact linear or nonlinear mode n for which the difference-frequency 
contact resonance occurs as n = p. The dominant difference-frequency component in 
Eqs.(l) and (1 1) is thus rj cp = r] cp diff = C cp ,diff and is given by Eq.(35) for n = p as 


C cp,diff GpCOs[(fi) c co s )t </) cc + fiss + Pcs $cs ^ • ( 50 ) 

where G p , given by Eq.(36), and the phase terms in Eq.(50) are obtained from Eqs.(15)- 
(18), (21)-(25), (37)-(43). It is important to point out in considering these equations that 
while the difference-frequency resonance frequency {co c -co s ) in RDF-AFUM is usually 
set to correspond to the lowest contact resonance mode of the cantilever (although a 
higher modal resonance could be used), the cantilever driving frequency oo c and 
ultrasonic frequency co s generally are set near (but not necessary equal to) higher contact 

resonance modes n = q and n = r , respectively, of the cantilever. For relatively small 
difference-frequencies, it may occur that q = r. Thus, the cantilever stiffness constant kc n 
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is appropriately given as kq, when involving the difference- frequency terms in Eqs.(15)- 
(18), (21)-(25), (37)-(43), the stiffness constant kcq when involving the cantilever drive 
frequency co c at or near the frequency of the qth cantilever resonance mode, and kcr when 

involving the ultrasonic frequency ro s at or near the frequency of the rth cantilever 
resonance mode. If co c and co s are not set at or near a contact resonance modal frequency 

of the cantilever, then it may be necessary to include more than one term in Eq.(l) and 
(11) corresponding to different values of q and r. 

It is seen from Eq.(36) that for a given value of {co c -co s ) the maximum value of 

Ccp,diff ideally occurs for a value of zo such that F"(zq ) is maximized. An examination 
of the force-separation curve of Fig.2 suggests that F"(zq ) is maximized near the bottom 
of the curve. It is important to note, however, that F'(z 0 ), while relatively small 
compared to that of the hard contact regime, is generally not equal to zero at that point. 
Strictly, the values of F"(zq ) and F'(zq) for a given zo are each dependent on the exact 
functional form of F(z 0 ) . A functional form for F(z 0 ) sufficiently quantitative to 
quantify F"(zq) and F'(z 0 ) is not typically available. However, experimental curves 
for F(z 0 ) can be obtained and compared 9 to the experimental curves of C C p,diff plotted 
as a function of zo. It is generally found that for a given difference-frequency (a> c -co s ) 
the maximum value of Ccp,diff occurs when F'(z 0 ) is negative, that is when z 0 is 
slightly smaller than the value of zo corresponding to the minimum value of F(z 0 ) in 
Fig.2. An examination of Eq.(36) suggests that a more exact approach to maximizing 
Ccp,diff would be not only to vary zq but also to vary slightly the difference-frequency 
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from the free space resonance condition until an optimal setting for both zo and the 
difference-frequency is achieved. 

The equations for G p and the phase terms in Eq.(50) may be simplified by using 
the values of the parameters estimated in Section III. All terms in Eqs.(15)-(18), (21)- 
(25), (37)-(43) involving the sample mass m s may be dropped to an excellent 
approximation. For ultrasonic wave and the cantilever drive frequencies in the low 
megahertz range we obtain, setting Aa = {co c -co s ) , that 


'cs 


tan 


-i rA Ae °) 


(51) 


, „ _i (r c k s + rs k cp)( Aa )-r s m c( Ao} ) +F '( z o)(r c + r s Ma) 

(f) cs « tan ~ Y 

Kp k s - OAv + YcYs )(Affl) + F (z 0 )\k +k s -m c (A co) ] 


(52) 


$cc ~ hffl 


-1 (Yc k s + Ys k cq H ~ Ys m c a b + F ' ( z 0 Xy c + Ys H 
^cq^s ~ (™ c k s YcYs) a> c F ( z oX^a/ ~ m c<°c ) 


(53) 


(f> ss ~ tan 


-1 (Yc^s Ys^cr^^s ys m c®s + Y s)® s 

kcr^-s ~ ( m c k s YcYs^s ^ ( z o)Acr ~ ) 


(54) 


and G p is given by Eq.(36) where 
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f 3 - * {ft 2 + r?(Affl) 2 } l/ 2 {[M« - (Affl)VA + r c r s ) 

^CS 

+ F(z 0 )(^ + *, - m c (A*>) 2 )] 2 (55) 

+ [(A<y)(/Ap + /A) - (Ao) 3 r 5 m c + F(z a )to c (r s + y c )] 2 V l/2 , 

Qcc * ^ c { A + ^ A )] 2 + rk 2 } 1 / 2 {[V, 

— co c (m c k s + y c y $) + F (z 0 )(k C q + k s — tn c Q) c )] (56) 

K (yA C q + Yaks') ~ ®cYs™c ^ AH' A /c)] } > 

a, * F,{[*, + F(z 0 )] 2 + 7> s 2 } 1/2 {[U 

— <y 5 (m c k s + y c y s ) + F (z 0 ){k cr + k s — tn c o ) s )] (57) 

s^cr Yc^s) ~ Y s m c k' (Zo)&c(.Ys 7c)] } » 

2cs' * ( z o ) { I k cr k s — co s (m c k s + YcYs) ( z o)(k C r ^ k s — tn c co s )] 

(58) 

\-*®s(.Yskcr Yc^s) ~ ®sYs™c {z 0 ')a> s (y s + y c )] } , 

and 

Qsc * ~~ k" c F' (z Q ) { \k C qk s — co c (m c k s + y c y s ) + F ' (z Q )(k C q + k s — tti c co c )] 


( 59 ) 
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\®c (.Ys^cq Yc^s) 0} cYs m c F* (Zo)(Dc(.Ys /c)] } 

The phase term T in Eq.(50) is given by Eq.(43) and is quite complicated. 
However, advantage can be taken of the fact that k s is generally quite large compared to 
other terms in the numerators of Q cc , Q ss , Q cs , and Q sc ; the denominators of these terms 
are very roughly all equal. Hence, the magnitudes of Q cc and Q ss are usually quite large 
compared to those of Q cs and Q sc . The terms involving the product Q cc Qss thus dominate 
in Eq.(43) and we may approximate T as 


i a cc - a ss = tan 


Ys m c 


K + 


- tan 


-1 


Yc a s 


^cr + (, z o) m c a s 


(60) 


where a cc and a ss are obtained from Eqs.(24) and (25), respectively. To the same extent 
that T may be approximated by Eq.(60) we may approximate G p as 


„ F"(z 0 ) D cs 

~ n r. zZCCZZSS 


R 


cs 


(61) 


It is seen from Eqs.(50)-(61) that both the amplitude and phase of the difference- 
frequency signal C cp, diff are dependent on F s , F c , k s , k c , y s , and y c in addition to (O c and 

co s . Since k s is proportional to the Young modulus of the material, the dependence of 
Ccp,diff on 7s an< i k s means that scans of the sample contain information about the elastic 


stiffness of the sample as well as information about surface damping, hence viscoelastic 
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properties of the sample surface. Subsurface features of the sample are obtained via the 
dependence of the difference-frequency signal amplitude on F s and via the dependence of 
the difference-frequency phase signal on 0, since both F s and 0 vary as the result of 
ultrasonic wave scattering from subsurface features. The signal response for F1FM is 
generally given by the same equations as those for RDF-AFUM except that a single mode 
p may not necessarily dominate the signal, if the difference-frequency is above the lowest 
contact resonance frequency of the cantilever. A sum of the largest modal contributions 
is thus calculated for HFM to obtain the signal output. However, the difference- 
frequency in HFM generally is set well below the lowest contact modal frequency of the 
cantilever. In this case the appropriate equations are identical to those of RDF-AFUM 
with p equal to the lowest linear or nonlinear contact modal frequency of the cantilever. 

B. Ultrasonic force microscopy 

In ultrasonic force microscopy (UFM) the cantilever drive frequency co c and drive 
amplitude F c are zero; the surface drive amplitude F s and the drive frequency co s of the 
wave generated by the transducer at the bottom of the sample are nonzero. UFM can be 
operated at quite large frequencies, even in the gigahertz range. Although the vibrational 
response of the cantilever is certainly quite small at such frequencies, operation at a tip- 
surface separation distance z 0 corresponding to the nonlinear regime of the force- 
separation curve, where F"(z 0 ) is maximum, will produce a detectable static or “dc” 
signal from the interaction nonlinearity. The generated static signal is called the 


ultrasonic force. 
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The nonlinear force-separation interaction results in a static displacement of the 
cantilever rj c stat given as 


Vc^tat S Vcn^tat 


(62) 


where rj cnstat is the contribution from mode n given by 


Vcn^tat £ cn + Ccn,stat 


(63) 


and e cn and Ccn,stat are given by Eqs.(13) and (20), respectively. Terms in Eq.(20) 

involving Q cc and Q sc are zero, since F c is zero for UFM. We assume operation of the 
UFM in the nonlinear regime where F"(z 0 ) is maximized andF'(z 0 )is negative. Using 
the values of the parameters estimated in Section III in the megahertz range of 
frequencies, we may approximate the nonzero terms Q ss and Q cs in Eq.(20) by Eqs.(57) 

and (58), where k cq is replaced with k cn . We obtain 


tfcn.stat 


K 

kcn^s “*■ F ( z o)(^cn + ^s) 


{F{z 0 ) 


+ ^“^[2 £q + Qcs + Qss + 2 QcsQss C °S« 55 ]} 


(64) 


4 
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where Sq is given by Eq.(21) and a ss is given by Eq.(25). To the extent that Q ss is 
much larger than Q cs because of the occurrence of k s and y s co s in the numerator of 
Q ss , Eq.(64) may be simplified by dropping the terms involving Q cs . 

Eq.(62) admits all cantilever modes as contributors to the magnitude of the UFM 
signal. However, Eq.(64) shows that the contribution to T] c stat for a given mode n is 

dependent on k cn such that for both the ultrasonic and non-oscillatory contributions to 
Eq.(64) an increase in k cn results in a decrease in the magnitude of the contributions for 
that mode. Since k cn increases in magnitude with increasing n, the contribution to rj c stat 

from a given mode generally decreases with increasing mode number for both the 
ultrasonic and non-oscillatory components of ri cn stat , although the exact relationship is 

highly dependent on the values of y c , y s , k s , m c , and co s that appear in Eq. (64). 

The dominant contributions from the second term on the right-hand side of 
Eq.(64) for a given ultrasonic drive frequency co s occur for those cantilever modes 

having values of \ k cn + F'(zq)] near the value m c co s . The largest contributions occur for 
values of co s near a contact modal resonance frequency of the cantilever 
yl\k cn + F'(zq)\I m c . In contrast, the first term on the right-hand side of Eq.(64) and the 
component of the second term involving f 0 are independent of frequency and thus make 
the major contributions when the ultrasonic drive frequencies are in the gigahertz range. 
These terms predict that a static signal exists even without the presence of an ultrasonic 
wave propagating through the sample and results directly from the interaction of the 
cantilever with the sample surface via the interaction force, as would be expected. 
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It is seen from Eq.(64) that rj cn stat is dependent on both F s and k s . This means 

that scans of the sample contain information about the elastic stiffness of the sample 
through k s as well as information about subsurface features via the dependence of the 
amplitude on F s . The dependence on y s means that UFM is sensitive to the viscous 
properties at the sample surface as well. 

C. Atomic force acoustic microscopy and force modulation microscopy 

Both for atomic force acoustic microscopy (AFAM) and force modulation 
microscopy (FMM) the cantilever drive amplitude and frequency are zero. As in UFM, 
the surface drive amplitude and frequency <n s are nonzero. However, unlike UFM, the 
surface drive frequency is limited to a range of frequencies that produces measurable 
displacement amplitudes of cantilever oscillation. In contrast to UFM, the tip-surface 
interaction distance is set to operate in hard contact, the “linear detection regime” of 
operation, where z 0 is small, F(z 0 ) is repulsive, F'(z 0 ) is large and negative, and 

F"(z 0 ) is negligible. In the “linear detection regime” no difference-frequency or 
harmonically generated signal is detectable, since for hard contact F"(zq) is effectively 
zero. The cantilever displacement amplitude rj cn i in corresponding to the nth mode is 
then obtained from Eq.(45) as T/ cn j in = e cn + Ccnjin ■ Th £ contribution C cn ,li n resulting 
from the nonlinearity is given by Eq.(26) and is seen to be zero, since F”(z G ) is 
effectively zero in the “linear detection regime.” The remaining contribution e cn to the 


cantilever displacement amplitude is given by Eq.(14). For values of the parameters 
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estimated in Section III we may approximate e cn , hence t] C n,lin ■> i n the low megahertz 
range of frequencies as 


Wcn,lin ~ ^cn * Qcs COS (fO s t (j) ss + 0) (65) 

where Q cs is given by Eq.(58), <f> ss by Eq.(54), and 9 .by Eq.(4). 

Note that both Q cs and (j) ss depend on the magnitude of F'(z 0 ) . For sufficiently 
hard contact F'(z 0 ) becomes very large and negative and may dominate the terms in 
Eq.(65). Under such conditions we obtain 


and 


Qcs ~ F s {(k cn + k s nifXOg) + {/ c + G > s } 


$ss 


tan 


(r c + r s )<»s 

ken k s ~ m c (O s 


( 66 ) 


(67) 


Eqs.(65)-(67) show that both the amplitude and phase of the cantilever 
oscillations depend on kc n , k s , y c , y Sj and co s . For AFAM driving frequencies oo s near a 

cantilever resonance «Jk cn / m c corresponding to mode n, the signal amplitude is large 

and the cantilever displacement is dominated by that mode. Although Eqs.(66) and (67) 
appear to be independent of F'(z 0 ) , the equations are subtly dependent on F'(z 0 ) via the 

restrictions F'(zq) places on the resonance modes available. As pointed out in Section 
II.B.zh, larges negative values of F'(z 0 ) eliminate from consideration all modes for 
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which F'(zq) < -k cn . For FMM, oo s is much smaller than ^jk c \ / m c , the fundamental 

cantilever resonance frequency, although the fundamental resonance may not appear in 
the calculations if F'(zq) < k c \ . 

D. Intermittent contact mode 

The intermittent contact mode (TappingMode) is a standard feature on many 
atomic force microscopes in which the cantilever is driven in oscillation, but no surface 
oscillations resulting from bulk ultrasonic waves are generated (i.e., F s and co s are zero). 
Thus, the intermittent contact mode cannot be used to image subsurface features, but 
interesting surface properties and features can be imaged. Since intermittent contact 
mode can be used in both hard and soft tip-surface contact (i.e. the linear and nonlinear 
regimes, respectively, of the force-separation curve), the cantilever displacement rj cn j in 

for mode n is given most generally as 

1 lcn.lin = %cn + £ cn,lin (68) 

where £ cn is given by Eq.(14) with the term involving Q cs set equal to zero and C cn ,li n is 
given by Eq.(26) with all terms involving Q cs and Q ss set equal to zero. 

For the soft contact regime the expression for rj cn i in is 

Vcn,lin = H cos (0) c t - (j} cc + A) (69) 


where 
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A = tan -1 sin (& + Fee <t>cc a cc) 5 (70) 

cos (J3 C + /j. cc — (j> cc — cx cc ) + ( Qcc 

W = -^e 0 F''(zo)(Ql + Ql + 2Q cc Qs C cosa cc ) v2 , (71) 

^ cc 

and 

H = [Ql + W 2 + 2 Q cc Wcos(/3 c + Mcc -</> cc - a cc )f 2 (72) 

where Q cc is given by Eq.(56), Q sc by Eq.(59), ^ cc by Eq.(53), n cc by Eq.( 27), s 0 by 
Eq.(21); a cc , fd c , D c , and R cc , are given by Eqs.(24), (29), (31), and (34), respectively, 
with the terms involving m s set equal to zero. 

The complexity of the cantilever response Tj cn j in is greatly reduced for the hard 
contact regime, where F"( z o) is negligibly small and F'(zq) is very large and negative. 
For hard contact A and a cc are approximately zero and we obtain 


Gendin ® Qcc COS {C0 c t $ C c) 


(73) 


where 


99 99 — 1/9 

Qcc ~ FcKKn k s ~ Vc ) Oc 7s) ®c ] 


(74) 


and 


0 CC = tan -1 + . 

ken + k s — m c (O c 


( 75 ) 
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The dependence of i) cn j in on the material damping coefficient y s and the sample 

stiffness constant k s , both for the hard and the soft contact regimes, means that the 

intermittent contact mode can be used to assess the viscoelastic properties of the material 
irrespective of the regime of operation. 

V. IMAGE CONTRAST 

All the above equations, except for Eqs.(4) and (6), were derived for constant 
values of the cantilever and material parameters. If, in an area scan of the sample, the 
parameters remain constant from point to point, the image generated from the scan would 
be flat and featureless. We consider here that the sample stiffness constant k s may vary 

from point to point on the sample surface. Since k s is proportional to the Young 
modulus E , this means that E also varies from point to point. We assume that the value 
of the sample stiffness constant k' s at a given point on the surface differs from the value 

k s at another position as k' s = k s + Ak s . For any function f(k s ) having a functional 
dependence on k s , a variation in k s generates a variation in f(k s ) given as 

Af = ( df /dk s ) 0 Ak s , where the subscripted zero indicates evaluation at k s . A similar 
expression can be obtained for the material damping parameter y s , but we shall not 
consider such variations here. 

A variation in k s produces a variation both in the amplitude and phase of the 

signal generated by the cantilever tip-sample surface interactions. The variations in 
amplitude and phase can be used to generate amplitude and phase images, respectively, in 
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a surface scan of the sample. We first consider images generated by the phase variations 
in the signal. 

A. Phase-generated images 

The phase factors involved in RDF-AFUM and HFM are given from Eq.(50), (4), 
and (5) to be (j) cc , <j> ss , /3 CS , <j> cs , T , and % ; the phase factors involved in AFAM and 
FMM are, from Eq.(65), <f) ss and % ; the phase factors involved in the intermittent contact 
mode are, from Eq.(69), (f > cc , and A . Each of these phase factors is dependent on k s and 
the variations in the phase factors resulting from variations in k s are responsible for 
image generation when using phase detection of the A-AFM signal. The exact 
dependence of the phase on k s , however, is different for soft and hard contact regimes. 


i. Soft contact regime 

For the soft contact regime the appropriate variations in the phase factors relevant 
to F1FM and RDF-AFUM are 


A P cs = 


r d(3 cs A 

V dk s jq 


Ak s = - 


r s Aa 


[k s + F'(z 0 )] 2 + rl(Aco) 


2 


(76) 


where 


A <i>cc 


^cc 

^cc 


A k s 


(77) 


dec — irsKq 7 F 1 ( z o )/ sk C q + F' ( z 0 ) (7c 7s)]®c 


(78) 
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[7c Ys 2 7s m c ( k cq F (zq)]o) c + m c y s (O c 


and 


where 


where 


and 


3,2 


Bcc {[/c^s Ys k cq Y s)\®c Ys m c G) c } 


( 79 ) 


^ n*o)(*« Vc) YcYs^c ) > 


A <fi ss =-^-Ak s 


B 


( 80 ) 




A ss = [Ys k cr + 2F, ( z o)7s k cr + ^'OofOc + / S )K 


( 81 ) 


and 


+ - 2 7s m c( k cr + F ( z 0 )]®S + 


3,2 


{[/c^y Ys k cr F (zq )(/ c 4- ^s-Xk^y Ys™c®s } 


( 82 ) 


and 


+ {[ k cr m c®s F ( z o)] k s F {. z o){.k cr m c CO s ) YcYs^s } > 


Wcs=~i^M s 

**cs 


( 83 ) 


4s = tMcp + 2F '( z o)7s k cp + F ( z o) 2 (Yc + ^)K A ®) 


( 84 ) 


+ 1 Jc 7s ~ 2 7s m c ( k cp + F ( z o)]( A ®) 3 + ™cYs( A a>) 5 


B cs {\.7c k s + 7s k cp + F ( z o)(/c + 7s)]( Aa> ) 7s^c ( A< ^) } 


( 85 ) 
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+ {[k C p - ™ c ( A(i >) 2 + ^OoXfc: + ^Oo )\- k cp - w c (Ao) 2 ] - / c / 5 (A®) 2 } 2 . 


To the extent that T = a cc -a ss , as given by Eq.(60), we may write 


A T = A a cc 


Ys^c 

ft + F'(z 0 )] 2 +y s 2 ffl| 


Ak s 


( 86 ) 


The phase term Ax is given by Eqs.(6) and (7). 

The appropriate variations in the phase factors relevant to the intermittent soft 
contact mode are A a cc , A<f > cc , and AA . The factor AA is obtained from Eq.(70) as 


AA = 


1 (Qcc / ^0 COS(/7 c + fj-cc &CC &CC ) 

[cosO? c + n cc - (/>cc -a cc ) + (Q cc / W)\ + sin (J3 C + H cc — (free ~ ^cc ) 


x (A(3 C + Afj. cc - A (free - Aa cc) 


(87) 


where 


*0C=- r^2 . ( 88 ) 

Ys®c 

A(f> cc \s given by Eq.(77), and Afi cc is obtained from Eq.(27). To the extend that Q sc is 
much smaller than Q cc , we get from Eq.(27) that A/u cc = Aa cc where A a cc is given by 


Eq.(86). 
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ii. Hard contact regime 

For the hard contact regime where F'(z 0 )is very large and negative, the relevant 
phase variations are obtained from Eqs.(67) and (75) as 


(k s + k cq m c m c ) + (/c + /,$■) 


Afe= . 

(*, + k cr m c a > s ) + (y c + y s } co s 


Eqs.(89) and (90) are appropriate to AFAM and FMM modalities as well as to the 
intermittent hard contact mode of A-AFM operation. 

iii. Dependence on the Young modulus 

Hertzian contact theory provides that the sample stiffness constant k s is related to 
the Young modulus E of the sample as 14 




(91) 


where v is the Poisson ratio of the sample material, E T and v T are the Young modulus 
and Poisson ratio, respectively, of the cantilever tip, and r c is the cantilever tip-sample 


surface contact radius. Hence, 
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Eq.(92) can be used with Eqs(76)-(90) to ascertain the fractional variation in the Young 
modulus AE / E from measurements of the phase variation in the signal from an 
appropriate A-AFM modality. 

B. Amplitude-generated images 

The amplitude G p of the RDF-AFUM signal is given by Eq.(61) to a good 

approximation for most applications. The fractional variation in the signal amplitude 
A Gp/Gp resulting from variations in the sample spring constant k s , hence Young 

modulus E , makes a considerable contribution to the image contrast when operating in 
the amplitude detection modality. The fractional variation in amplitude is 

^/L = _L (??£_) Ak 

G p G p { dk s ) 0 

(93) 

t i (e&,) , i (BQcA il M 

|(6ss v dk s Q cc v 8k s D cs dk s v R cs v Jq|j 


i rga s i o) 

QssVdK ) 0 [(t s +F(z 0 )) 2 + i ] 1/2 


where 
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( 94 ) 

( a ssk s ^ss)^ss (c ss k s + d ss )c ss 

( a ssk s k*ss ) ( c ss k s ^ss ) 


a ss = k cr ~ m c<°s + F ( z 0) ’ ( 95 ) 

b ss — F {zQ^){k cr — m c co s ) — Y cY s®s > ( 9 ^) 

^ss ~ Y c®s ’ ( 9 < 7 ) 

^ss ~ Yskcr a s ~ m cYs ( °s ^( z o)(/c Fv) > ( 9 ^) 


i pgcc] ^+n z o) 

Gccl^Jo [(^ + F'(z 0 )) 2 + C c 2 c ] 1/2 

(99) 

_ (®CC k S k>cc)^CC (pcc k s decree 

(decks ^cc ) (pecks d C c ) 

«CC = k cq - m c a>e + F'(z 0 ) , (100) 


k cc F'(zo)(k C q m c O ) c ) YcYs®c ’ 


(101) 
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^cc Yc a c 


( 102 ) 


d cc Ys^cq^c ™cYs®c F ( z o)(Yc Ys) 


(103) 


and 


Res 

1 

CD 

( D ^ 

^ CS 

1 ** 

D cs 

dk s \ 

\Rcs J 

Jo [ft 2 +4 f 1 


(a cs k s + b cs )a ss + ( c cs k s + d cs )c cs 

( a cs k s + &cs ) i^cs^s ^cs) 


(104) 


^CS k cp m c (A©) + ^ ( Z 0) ’ 


(105) 


= F'(z 0 )(k cp - m c (Ao)) 2 ) - y c y s (Aco) 2 , (106) 


d cs Yc (Aft?) j 


(107) 


*cs = Ys k cp( A(0 ) - Wcr s (A«) 3 + F(zq)(j c + r 5 ) ■ (108) 


It is apparent from Eqs.(93)-(108) that, although the RDF-AFUM signal 
amplitude per se is highly dependent on F r '(z 0 ) and on the cantilever and ultrasonic 


drive amplitudes F c and F s , respectively, the magnitude of the fractional variation 
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A G p /G p in the RDF-AFUM signal amplitude resulting from variations in the sample 
spring constant k s is independent of F c , F s , and F"(zq). However, A G p /G p is 
dependent upon the values of the cantilever spring constant k cn where n = p, q, and r as 
discussed in Section IV. A. The values of k cn in turn are highly dependent on the choice 
of cantilever and the frequency chosen to drive the cantilever into resonance. Although 
A G p /G p makes a considerable contribution to image contrast, it is not the only 

contribution. As with all A-AFM techniques the resolution of the image digitizer, the 
dynamic range and signal-to-noise features of the electronic components, the sharpness of 
the cantilever tip, and the bonding of the ultrasonic transducer among other factors also 
contribute to the image contrast. The contrast for RDF-AFUM, however, cannot 
generally exceed that rendered by A G p /G p . The magnitude of the signal variation for 

HFM is given by the same equation as for RDF-AFUM except that a single mode may 
not necessary dominate the signal. A sum of the largest modal contributions may be 
appropriate to calculate for HFM the cantilever displacement. 

The amplitude of the UFM signal is given by Eqs.(62)-(64). Assuming that Q cs is 

small compared to Q ss and Sq, hence negligible in the calculations, we obtain the 
fractional variation in the cantilever displacement amplitude for the nth mode to be 


Vcn^tat 




cn 




1 


^s^cn^s + FitoWcn ^y)l F(zq) + (2sq + Q$ s 


} Ak s 

'Jo 


( 109 ) 
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where Q ss is given by Eq.(57) with q = n, (8Q ss /dk s ) oby Eq.(94) with q = n, by 
Eq.(21), and (5f 0 /5^) 0 by 


= F( z o)kcn 

v dk s ) 0 k cn k s + F(z 0 )(k cn + k s ) 


( 110 ) 


For AFAM and FMM the cantilever displacement amplitude is from Eq.(65) 
dependent on Q cs where Q cs is given by Eq. (58). The fractional change in the signal 
amplitude for mode n is obtained to be 


^ Qcs 
Qcs 


ken k s m c (Q s 


(k C n k$ m c a) s ) + (/ c + Y s ) 


2 2 


(111) 


The availability and dominance of modes are discussed in Section IY.C. 

For the intermittent hard contact modality the amplitude is dependent on Q cc 
which for hard contact is given by Eq.(74). The fractional change in the amplitude for a 
given mode n is obtained as 

^ Qcc k cn + k s — m c 6) c a 7 

/_ /y 9 9 

Qcc (k cn + k s — m c o ) c ) + (y c + y s ) co c 


(112) 
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It is interesting to note that the variation in amplitude for intermittent hard contact is 
identical to that of AFAM and FMM except that in intermittent hard contact the drive 
frequency is co c whereas in AFAM and FMM the drive frequency is co s . 


VI. APPLICATION TO LARC-CP2 POLYIMIDE FILM 

The above equations are applied to the assessment of AE/E in LaRC-cp2 
polyimide polymer from measurements 9 of the variations in the phase signal using RDF- 
AFUM. The values of the relevant material and cantilever parameters are 9 k s = 96.1 N 
m' 1 , k cl = 14 N m' 1 , y s = 4.8 x 10' 5 kg s' 1 , m c = 3.9 x 10' 12 kg, E = 2.4 GPa, a/2 = 12.7 
pm, a = 85 m" 1 , F'{zq)= - 53, co c Hn = 2.1 MHz, co s l2n = 1.8 MHz, and Acolln 

0.3 MHz. From Eq.(50) the variation in the phase signal is given as 
- A(j) cc + A(3 CS - Afl cs + AT + Aj) . From Eqs.(6),(7),(76)-(86) we calculate that 
A(3 CS and AT = Aa cc make by far the dominate contributions to the RDF-AFUM phase 
signal variations. The contribution from Ay is relatively small because the thickness of 
the sample a/2 does not correspond to a resonance thickness for the ultrasonic wave. The 
factors A<j> ss and A<j) cc are relatively close in magnitude and result in relatively little net 
contribution to the phase variations in Eq.(50). The factor A <j> cs is also calculated to be 
small with a sign opposite to that of (A<f> ss - A<f> cc ) and results in further minimizing the 
net phase contribution from these three terms. 

From the above equations and a measured phase variation of 13.2 degrees 
obtained in the RDF-AFUM phase image 9 we calculate a value of approximately 24 
percent for the variation in the Young modulus for the material. This value is also in 
good agreement with a value of roughly 21 percent obtained from independent 
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mechanical stretching experiments in which the increase in the modulus is attributed to 
the growth during stretching of a crystalline phase having a larger Young modulus than 
that of the original amorphous phase 20 . The growth of the crystalline phase in LaRC-CP2 
is attributed to internal stresses generated during curing by the different thermal 
expansion coefficients of the polymer material and gold nanoparticles embedded in the 
polymer matrix. 

A further test of the present model can be obtained from the intermittent soft 
contact image taken concurrently with the RDF-AFUM image of the LaRC-CP2 
specimen 9 . The variation in the phase signal is measured from the intermittent soft 
contact micrograph to be roughly 1.5 degrees. From Eq.(69) the variation in the phase 
signal for intermittent soft contact is analytically given as (-A^ CC + AA). Using the 

above-stated values of the material and cantilever parameters in Eqs.(86)-(88) ,we 
calculate that for the measured variation of 1.5 degrees in the intermittent soft contact 
image the variation in the Young modulus is roughly 18 percent. This value is in good 
agreement with the values obtained from the RDF-AFUM image and from the 
independent, mechanical stretching measurements. 

VII. CONCLUSION 

The various modalities of acoustic-atomic force microscopy (A-AFM) have 
become important nanoscale characterization tools for the development of novel 
materials and devices. Most of the information obtained from A-AFM has been 
qualitative because of the lack of a comprehensive analytical model to render the data 
quantitative. The most significant impediment to the development of such a model has 
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been the nonlinearity of the cantilever tip-sample surface interaction force. We have 
developed a detailed mathematical model of this interaction by assuming that the 
interaction is appropriately represented by a nonlinear spring coupled at the opposite ends 
by linear springs representing simple harmonic oscillators for the cantilever and sample 
surface. The dynamics of the coupled springs are described by a pair of coupled 
differential equations that are solved using a standard iteration procedure. Only flexural 
vibrations of the cantilever and out-of-plane oscillations of the sample surface are 
considered in the present derivation. 

Solutions are obtained for specific A-AFM modalities including the commonly 
used intermittent contact mode (TappingMode), force modulation microscopy (FMM), 
atomic force acoustic microscopy (AFAM), ultrasonic force microscopy (UFM), 
heterodyne force microscopy (HFM), and resonant difference-frequency atomic force 
ultrasonic microscopy (RDF-AFUM). Image generation and contrast equations are 
obtained for each of the aforementioned A-AFM modalities assuming for expediency that 
the contrast results only from variations in the sample stiffness constant. Since the 
sample stiffness constant is related directly to the Young modulus of the sample, the 
contrast can be expressed in terms of the variation in the Young modulus from point to 
point as the sample is scanned. 

A portion of the present solution set was obtained previously 9 using feasibility 
arguments to render the analytical equations and applied to the assessment of the 
fractional variation in the Young modulus AE / E of LaRC™-CP2 polyimide polymer 
from measurements of the variations in the phase signal using RDF-AFUM. The present 
derivation provides a solid analytical framework for the feasibility arguments. The 
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calculation of 24 percent for A El E is in good agreement with a value of 21 percent 
obtained from independent mechanical stretching experiments in which the increase in 
the modulus is attributed to the growth during stretching of a crystalline phase having a 
larger Young modulus than that of the original amorphous phase 20 . 

Further application of the model to intermittent soft contact images of the same 
polymer material yields a value of 18 percent for the variation of the Young modulus. 
This value is also in good agreement with the mechanical stretching experiments. The 
difference between the RDF-AFUM and intermittent soft contact assessments of AE/E is 
attributed to error propagation resulting from the relative complexity of terms, the 
approximations used, and the number the terms encountered in the respective model 
calculations. Nonetheless, the agreement between the RDF-AFUM and the intermittent 
soft contact assessments is quite good considering that the calculations are performed on 
rather disparate equations. More importantly, the agreement between the model 
predictions and the measurements from independent mechanical tests provides quite 
strong evidence for the general validity of the present model. 

The present model can also be used to quantify the image contrast from variations 
in the sample absorption coefficient y s or from a combination of absorption coefficient 

and Young modulus variations in the material. Space limitations prohibit the inclusion of 
such contrast mechanisms here, but the effects can be derived straightforwardly by the 
reader from the equations derived above. Although the present model is developed for 
flexural oscillations of the cantilever and out-of-plane vibrations of the sample surface, 
the model can in principle be extended to include other modes of cantilever oscillation 
and sample surface response. It is anticipated that such a development would provide 
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even greater opportunities for obtaining quantitative information on material properties 
using the various A-AFM modalities. 
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FIGURE CAPTIONS 

Fig. 1 . Schematic of cantilever tip-sample surface interaction: Left side shows relative 

positions of the cantilever tip and sample surface in the presence of an interaction 
force. Right side shows spring model representing the dynamics of the tip-surface 
interaction, zo is the quiescent tip-surface separation distance, z the oscillating 
tip-surface separation distance, r/ c the displacement (positive down) of the 

cantilever tip , tj s the displacement of the sample surface (positive up), k cn is the 
nth mode cantilever spring constant, m c the effective cantilever spring mass, k s 
the sample spring constant, m s the effective sample spring mass, and F'(z 0 ) and 
F"(zq) are the linear and first-order nonlinear sample stiffness constants, 
respectively. 

Fig.2. Schematic of interaction force as a fimction of the separation distance between 
cantilever tip and sample surface. 
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Fig. 1. 
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QuickTime™ and a 
TIFF (LZW) decompressor 
are needed to see this picture. 


Fig. 2. 



